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基于 空气 网 络 法 的 管 路 系统 
及 阀门 计算 模型 
孙 科 ， 刘 高 文 ， 马 春田 *， 和 牛 嘉 嘉 ， 
(1、 西 北 工业 大 学 动力 与 能 源 学 院 ， 陕 西 西安 ，710072; 
2、 东方 汽轮机 有 限 公司 四 川 德阳 ，618000) 
摘 要: 空气 网 络 法 因 其 快速 实用 的 特点 ， 计 算得 到 空气 网 络 中 各 截面 各 种 参数 ， 在 航空 发 动机 空气 系统 的 计算 研究 中 得 到 
广泛 应 用 。 将 空气 网 络 法 的 思想 应 用 到 压缩 空气 供应 系统 中 ， 其 中 主要 有 了 阀门、 管道 、 局 部 损失 等 主要 元 件 。 因 原 有 程序 中 没 
有 阀门 计算 方法 ， 首 先 建立 阀门 的 计算 模型 ， 通 过 实验 得 到 阀门 的 流量 特性 ， 形 成 计算 方法 ， 写 入 空气 网 络 程序 中 。 设 计 了 一 
种 新 型 数据 结构 来 表达 空气 网 络 ， 重新 编写 了 面向 对 象 语言 程序 。 使 用 空气 网 络 程序 计算 不 同 阀门 开 度 工 况 下 管 路 的 流量 及 
要 截面 上 的 参数 ， 并 与 实验 值 进 行 对 比 ， 偏 差 很 小 ， 从 而 验证 计算 方法 的 准确 可 人 和信。 由 此 得 到 压缩 气体 供应 系统 的 空气 网 络 模 
型 ， 应 用 于 管 路 系统 调控 中 。 进 行 空 气 网 络 计 算 ， 根 据 实验 台所 需 的 气体 流量 及 状态 参数 得 到 所 需 阀门 的 开 度 ， 快 速 调 整 阀门 
到 所 需 位置 。 
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Abstract: Because of its fast and practical characteristics, the 1D network model has been widely used in the 
research of secondary air System in engine. The 1D network model is applied to the compressed air Supply System， 
which mainly has the main elements such as valves, tubes, losses and so on. Because there is no valve calculation 
method in the original program, the calculation model of the valve was established firstly, and the flow characteristic 
of the valve was obtained through the experiment, then the calculation method formed. A new data structure is 
designed to express the air network, and the object oriented programming language is recompiled. Using 1D network 
program to calculate the flow rate and the parameters of the main section of the system under different valve opening 
conditions, and the experimental results are compared with the computational results, the deviation is very small, so 
as to verify the accuracy of the calculation method. Calculation model of compressed air supply system which is 
applied to control the supply system. Get the required opening of valve which control system to provide suitable 
compressed air for experiment. So fastly adjust opening in experiment. 
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0 引言 


压缩 空气 是 航空 发 动机 工作 中 必 不 可 少 的 工 质 ， 
同时 也 是 模拟 实验 中 的 工 质 。 实 验 室 中 ， 压 缩 空气 供 
应 系统 为 实验 台 供给 所 需 质 量 、 品 质 的 空气 。 空 气 供 
应 系统 的 设计 ， 应 该 减 小 气流 损失 ， 并 且 准 确 无 误 为 
实验 台 服 务 。 而 阀门 作为 管 路 系统 中 调节 气流 的 关键 
元 件 ， 有 十 分 重要 的 作用 。 调 节 阀 门 的 开 度 ， 得 到 气 
流 在 系统 中 的 不 同 状态 。 通 过 空气 网 络 计 算得 到 阀门 
需要 打开 的 开 度 ， 从 而 迅速 调整 到 位 。 因 此 压缩 气体 
供应 系统 的 空气 网 络 计算 较 有 意义 。 国 内 外 学 者 应 用 
空气 网 络 法 进行 了 很 多 的 研究 ， 大 部 分 研究 是 将 网 络 
法 应 用 于 发 动机 空气 系统 中 。Yannickt 等 通过 开源 软 
件 CalculiX®@ 进 行 喷气 发 动机 二 次 空气 系统 模型 综合 
热机 械 分 析 ， 模 拟 二 次 空气 系统 中 流动 特性 及 热 结 构 
特性 的 综合 作用 效果 。 Benra 折 等 使 用 一 维 空气 网 络 模 
型 模拟 预 旋 系统 中 的 流动 状态 。 一 维 计算 的 结果 与 数 
值 计 算 结 果 和 实验 结果 进行 了 对 比 。 吴 丁 妆 局 提 出 了 
关于 内 流 系统 流动 换 热 的 网 络 计 算法 ， 计 算 稳定 性 和 
收敛 性 较 好 。 伐 升平 中 等 开发 了 用 于 模拟 发 动机 空气 
系统 非 稳 态 特 性 的 程序 。 郭 晓 杰 中 考虑 空气 系统 流体 
特性 (温度 、 压 力 与 流量 ) 和 固体 部 件 热力 特性 (温度 分 
布 ) 之 间 的 相互 作用 , 将 空气 系统 计算 与 固体 部 件 热 分 
析 结 合 起 来 ， 建 立 航空 发 动机 空气 系统 与 固体 部 件 2 
间 的 稳 态 流 - 热 耦 合 分 析 平 台 。 吕 亚 国 % 开 发 了 发 动机 
空气 系统 通用 分 析 软 件 ， 软 件 前 台 采 用 图 形 建 模 方式 
生成 网 络 计算 模型 ， 后 台 采 用 了 网 络 自动 识别 技术 。 
胡 肖 肖 中 等 借助 概率 思想 建立 相应 的 计算 连续 性 方程 
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和 能 量 方程 ， 建 立 相 应 的 随机 游 动 模型 ， 运 用 蒙特 卡 
罗 方 法 求解 各 节点 的 压力 ， 再 根据 流量 与 元 件 两 端 压 
力 的 关系 ， 计 算 流 经 各 元 件 的 流体 流量 ， 结 果 与 
Flowmaster 进行 了 对 比 ， 吻 合 较 好 。 潘 耘 峰 由 以 某 型 
航空 发 动机 为 例 ， 得 到 主要 元 件 的 计算 方法 ， 并 建立 
空气 网 络 模型 。 同 时 有 很 多 研究 ， 将 空气 网 络 法 应 用 
于 其 它 工作 中 。 陆 震 露 所 将 一 维 计算 方法 应 用 于 内 燃 
机 性 能 仿真 ， 并 开发 了 仿真 软件 。 张 光 鹏 、 许 诺 等 5 
应 用 Flowmaster 对 某 一 空调 系统 进行 了 数值 模拟 , 结 
果 表 明 Flowmaster 对 暖 通 空调 设计 有 着 积极 地 指导 意 
义 


本 文 将 空气 网 络 法 应 用 于 压缩 气体 供应 系统 中 ， 
建立 阀门 的 模型 ， 通 过 实验 得 到 特性 参数 ， 代 入 空气 
网 络 程序 中 。 进 行 空 气 网 络 计 算 ， 结 果 与 实验 结果 进 
行 对 比 。 通 过 空气 网 络 计算 可 快速 得 到 管道 中 重要 截 
面 的 气流 参数 。 

1 空气 网 络 及 计算 方法 
1.1 空气 网 络 模型 

根据 某 实 验 室 压缩 气体 供应 系统 中 管道 、 阀 门 的 
布置 ， 简 化 得 到 相应 的 空气 网 络 模型 。 图 1 为 某 实验 
室 压缩 气体 供应 系统 的 空气 网 络 模型 ， 该 空气 网 络 模 
型 中 有 一 个 进口 、 一 个 出 口 ， 六 个 腔 室 ， 十 个 分 支 。 
11 为 进口 (Inlet)， 了 为 出 口 (Outlet)。C1-C6 为 6 个 腔 
室 (Chamber)，E1-E43 为 元 件 (Element)。 网 络 中 有 管 
元 件 (Tube)、 局 部 损失 元 件 @Partloss)、 阀 门 元 件 (Value) 
这 三 种 , 共有 43 个 元 件 。 表 1 所 列 为 空气 网 络 中 各 个 
分 支 的 元 件数 量 和 分 支 的 进口 、 出 口 。 
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图 1. 一 维 空气 网 络 模型 


Fig 1. 1D network model 
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这 些 元 件 的 流 阻 和 传 热 特性 计算 模型 和 计算 方法 m 为 气体 质量 流量 ，2 为 气体 标准 状态 下 的 密度 ; 


的 建 并 可 借助 理论 推导 、 已 有 实验 数据 和 参考 文献 等 ， KK, 为 阀门 的 流量 系数 ， 定 义 是 温度 为 (5-40'C) 的 
并 通过 将 计算 结果 与 已 有 实验 数据 进行 对 比 来 对 计算 水 在 10Pa 压 降 下 ，1 小 时 内 流 过 阀 的 立方 米 数 ; 
模型 和 方法 进行 改进 和 完善 。 其 中 管 元 件 "、 局 部 损 7- 立 门 入 口 的 绝对 温度 ; 
失 元 件 中 的 计算 方法 在 已 有 的 空气 网 络 程序 中 已 G- 气 体 的 相对 密度 ; 
经 存在 ， 经 过 多 次 使 用 ， 已 经 得 到 了 验证 。 本 文 建立 亏 压 缩 系 数 ; 
阀门 元 件 的 计算 方法 ， 根 据 实 验 得 到 实验 室 所 采用 阀 大 绝热 指数 。 
门 的 流动 特性 ,之 后 计算 方法 写 入 到 空气 网 络 程序 中 ， 当 系 8 时 ， 非 阻塞 流 ， 进 口 总 压 计 算式 为 
进行 整个 空气 网 络 模型 的 工程 计算 。 
表 1 网 络 中 的 分 支 pi= ls Ze (1-b) 
Table1. Branchs in network model 人 
二 本 人 Number of Branchs Branch’s 式 中 yy 膨胀 系数 ， = 
elements inlet outlet XxX 7 
1 3 11 Cl 1.3 空气 网 络 数据 结构 
， x 大 型 空气 网 络 中 一 般 有 很 多 元 件 ， 并 且 有 多 个 进 
口 、 多 个 出 口 、 多 个 腔 室 、 多 个 分 文 ， 是 非常 复杂 的 
= 图 形 结构 ， 需 要 有 效 的 组 织 空气 网 络 数据 结构 才能 使 
4 2 Cl C2 计算 机 轻松 识别 。 文 献 [ 和 中 通过 定义 各 个 元 件 的 进出 
5 2 Cl C2 口 节 点 编号 ， 并 且 相 邻 的 两 个 元 件 间 ， 上 游 元 件 的 出 
, 3 口 节点 号 与 下 游 元 件 的 进口 号 下 烙 保持 一致 。 这 种 方 
法 的 原理 简单 ， 但 对 用 户 来 说 较为 繁 杀 ， 需 要 特别 熟 
悉 编号 规则 ， 并 且 当 网 络 中 减少 或 增加 一 个 元 件 时 ， 
8 9 C4 C5 几乎 所 有 的 节点 序号 和 元 件 序 号 都 需要 重新 定义 才能 
9 4 C5 C6 保证 分 文中 从 上 游 到 下 游 元 件 序 号 依次 增 大 的 顺序 。 
是 a 为 方便 使 用 提出 了 一 种 新 型 的 数据 结构 ， 如 图 2 中 简 
一 -一 易 的 空气 网 络 模型 ， 以 分 支 为 单元 ， 图 中 所 示 的 分 支 
1.2 元 件 计 算 方 法 间 的 有 向 连接 可 以 使 用 如 表 2 的 二 维 整形 数组 表示 ， 
数组 元 素 A[0][1]=1 表示 从 4 1 到 分 支 2 的 有 所 
闪 门 是 风 洞 系统 中 最 常用 的 空气 流量 调节 元 件 ， et i 


二。 ”实验 室内 门 调节 的 控制 参数 是 开 度 。 实 验 宝 中 所 采用 。 的 数据 结构 


电 ”的 阀门 是 电力 驱动 阀门 ， 单 座 套 简 调 节 阀 ， 流 开 型 。 、 
通过 阀门 的 流体 可 压缩 ， 因 此 进口 总 压 的 计算 存在 阻 /7 
塞 流 和 非 阻塞 流 两 种 情况 。 这 里 介绍 三 个 参数 ， Fk 为 让 了 全 
比 热 比 系数 ， 对 于 空气 取 1， 对 非 空气 介质 FH1.4: Ee 
陈 为 阀门 前 后 压 差 比 , 即 阀门 压 降 如 与 入 口 压力 py 由 SP 
于 管道 内 气体 马赫 数 很 小 ， 动 压 远 小 于 静 压 ， 因 此 静 -| 
压 近 似 等 于 总 压 ) 之 比 ; 各 为 临界 压 差 比 , 对 于 一 个 特 Va 
定 的 阀门 ， 当 产生 阻塞 流 时 ， 其 压 差 比 是 一 个 固定 党 Ma 
数 ， 对 单 座 套 简 型 流 开 形式 的 取 0.75。 NAN 
当 他 加 Fx 时 ， 出 现 阻塞 流 ， 进 口 总 压 计 算式 "9 AN 
为 图 2. 网 络 数据 结构 
0, TGZ Fig 2. Data structure of network 
D258K,Y rr (9 


式 中 O。 气体 标 准 状 态 下 的 体积 流量 ，QOs=mp， 


表 2 网 络 邻 接 矩 阵 
Table 2. Network adjacent matrix 


I 


AD 让 0 1 2 3 

J 
O(branch1) 0 0 0 0 
1(branch2) 1 0 0 0 
2(branch3) 1 0 0 0 
3(branch4) 0 0 0 


具体 到 每 个 分 支 的 数据 结构 ， 如 图 3 所 示 ， 每 个 
元 件 中 有 pnext 和 plast 两 个 指针 将 不 同 元 件 串联 起 
来 ， 组 成 双向 链表 。 通 过 面向 对 象 编程 中 类 继承 中 的 
思想 将 标准 类 派生 为 不 同 元 件 类 ， 从 而 不 同类 型 元 件 
能 够 组 成 一 个 线性 链表 。 只 需 按 照 顺 序 输入 各 个 元 件 
的 输入 参数 ， 无 需 对 节点 进行 编号 ， 计 算 机 自动 识别 
且 生成 链表 。 
Tube pnext 上 - 
hl plast 
图 3. 分 支 数据 结构 
Fig 3.Data structure of branch 
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具体 到 每 个 元 件 的 数据 结构 ， 如 图 4 所 示 ， 每 个 
元 件 类 中 有 函数 成 员 和 数据 成 员 ， 而 元 件 计算 函数 
Cal_ element(...) 从 标准 类 中 派生 而 来 ， 因 而 在 代码 编 


写 时 ， 不 同 元 件 计算 只 需 调 用 这 一 个 函数 即 可 。 另 外 
由 于 数据 成 员 也 定义 在 类 中 ， 在 程序 调试 过 程 中 方便 
监测 元 件 进 出 口 参数 的 变化 ， 而 原 有 数据 结构 中 ， 因 
为 不 同 气动 参数 定义 在 不 同 数组 中 ， 通 过 节点 序号 索 
引 该 元 件 参数 ， 在 调试 过 程 中 需要 反复 切换 数组 。 另 
外 链表 这 样 的 数据 结构 对 内 存 空间 的 节省 使 用 也 十 分 
有 利 。 


Cal element(...) 


图 4. 元 件 类 


Fig 4.Class of element 


2 实验 室 阀 门 流动 特性 实验 
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2.1 实验 条 件 、 测 量 技术 

实验 室 所 用 的 压缩 空气 来 自 螺杆 空气 压缩 机 ， 最 
大 供 气 压力 0.7MPa, 压缩 气体 在 气 钢 中 稳定 后 进入 管 
路 系统 。 
压力 测量 采用 压力 变 送 器 。 测 量 范 围 0-1MPa, 精 
度 等 级 0.2 级 。 传 感 器 显示 屏 输 出 压力 示 数 ， 另 外 传 
感 器 输出 4-20mA 电流 信和 号， 经 信号 隔离 器 转换 成 
1-5V 电压 信号 ， 电 压 信 和 号 由 研华 4711A 板 卡 采集 进 
入 工业 控制 计算 机 中 。 温 度 测 量 技术 介绍 ， 温 度 测量 
采用 K 型 热电 偶 ， 量 程 40 一 1350'C。 流 量 测量 技术 
介绍 ， 流 量 测量 使 用 标准 孔 板 节 流 装置 配合 差 压 变 送 
器 、 压 力 变 送 器 、 热 电 偶 组 成 。 差 压 便 送 加 精度 等 级 
0.1。 使 用 差 压 变 送 器 测量 得 到 和 孔 板 前 后 压 差 AP, 压力 
变 送 器 测 得 孔 板 前 压力 ， 热 电 偶 测量 得 到 气流 温度 。 
2.2 阀门 流动 特性 实验 

实验 室 采 用 的 电动 阀 为 单 座 套 简 阀 ， 流 开 形 式 。 
该 电动 阀 的 流量 特性 需 由 实验 得 到 。 实 验 中 ， 电 动 阀 
门 前 后 各 有 一 个 压力 传感器 , 在 阀门 上 游 处 布置 K 型 
热电 偶 ， 下 游 有 孔 板 流量 计 。 

由 式 (1-a) 得 ， 当 X>XrFk 时 ， 出 现 阻塞 流 ， 流 量 
系数 的 表达 式 为 


_ 0 cz 
” 2.58p rr 


二 


(2-a) 


由 式 (1-b) 得 ， 当 X<Xz 时 ， 非 阻 寒流， 流量 系 
数 的 表达 式 为 
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表 3. 量 纲 分 析 


Table 3. Dimensional analysis 


Symbol Result Symbol Result 
Qs LT k dimensionless 
Pi MTL’! XT dimensionless 
Ti 0 y dimensionless 
G dimensionless X dimensionless 
Z dimensionless 


流量 系数 表达 式 右 端 参 数 有 量 纲 ， 现 对 流量 系数 


的 量 纲 进 行 分 析 。 对 式 (2-a) 和 (2-b) 中 等 式 右 边 的 参数 
量 纲 分 析 结 果 如 表 3。 由 表 3 量 纲 分 析 结 果 及 式 (2-a) 


口 工 
门 = 


二 || 
[J 


和 (2-b) 分 析 得 流量 系数 的 量 纲 为 L173M -03 


阀门 不 同 开 度 下 ， 流 量 系 数值 不 同 。 实 验 中 ， 根 
据 阀门 前 后 压力 、 通 过 阀门 的 流量 、 阀 门 前 气流 温度 ， 
便 可 计算 得 到 阀门 不 同 开 度 下 ， 对 应 的 流量 系数 值 。 


3 空气 网 络 计算 结果 与 实验 对 比 


3.1 阀门 流动 系数 特性 曲线 

阀门 流动 特性 主要 随 开 度 变 化 ， 本 文 做 了 气 源 压 
力 0.7MPa、0.5MPa、0.35MPa 这 三 组 实验 。 如 图 5 
为 实验 结果 ， 图 中 3 条 曲线 分 别 为 3 个 不 同 气 源 压力 
时 的 流量 系数 曲线 。 三 条 曲线 相近 且 符 合 阀门 的 等 百 
分 比特 性 曲线 。 根 据 不 同 气 源 压力 下 各 开 度 对 应 流量 
系数 的 平均 值得 到 表达 式 y=157.480x*+45.915x+7.002 
代入 程序 计算 。 


air source pressure=0.7MPa 
一 一 一 airsource bressure=0.5MPa 
一 一 一 一 airsource pressure=0.35MPa 刀 
Sp 


0.2 0.4 0.6 0.8 1.0 
value 


图 5. 阀门 流量 系数 
Fig $5. Discharge coefficient of value 

3.2 网 络 计算 结果 与 实验 对 比 

对 管 路 的 空气 网 络 建 模 计 算 中 , 调整 阀门 的 开 度 ， 
形成 不 同 的 管 路 工 况 。 针 对 不 同 的 工 况 ， 编 号 空气 网 
络 程序 输入 文件 ， 并 进行 计算 。 
首先 是 气 源 压力 为 0.7MPa( 表 压力 ), branch5 中 电 
动 阀门 的 开 度 为 0.15 时 的 工 况 。 如 图 6 为 空气 网 络 程 
序 计算 结果 与 实验 测量 值 的 对 比 。 使 用 计算 与 实验 值 
相对 偏差 衡量 计算 准确 性 ， 相 对 偏差 的 定义 为 计算 值 
与 实验 值 差 的 绝对 值 与 实验 值 之 比 。 实 验 测 得 网 络 中 
四 个 点 的 压力 , 计算 与 实验 值 的 对 比 , 最 大 偏差 3.2%， 
最 小 偏差 0.2%。 计算 所 得 的 流量 值 与 实验 测量 结果 偏 
差 4.1%。 
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6. 沿 程 压力 ( 气 源 0.7MPa) 
Fig 6.Pressure along the way(0.7MPa) 
其 次 是 气 源 压力 为 0.7MPa( 表 压力 ), branch5 中 电 
动 阀门 的 开 度 为 0.8 时 的 工 况 。 如 图 7 为 空气 网 络 程 
序 计算 结果 与 实验 测量 值 的 对 比 。 在 阀门 处 的 压力 损 
失 比 开 度 为 0.15 时 小 , 计算 所 得 的 压力 值 与 实验 值 的 
对 比 ， 最 大 相对 偏差 0.6%， 最 小 相对 偏差 0.1%。 计 
算 所 得 的 流量 值 与 实验 测量 结果 相对 偏差 3.9%。 


1D network 


pressure(10Pa) 


图 7. 沿 程 压力 ( 气 源 0.7MPa) 

Fig 7.Pressure along the way(0.7MPa) 
其 次 是 气 源 压力 为 0.5MPa( 表 压力 ), branch5 中 电 
动 阀门 的 开 度 为 0.8 时 的 工 况 。 如 图 8 为 空气 网 络 程 
序 计算 结果 与 实验 测量 值 的 对 比 。 计 算 所 得 的 压力 值 
与 实验 值 的 对 比 ， 最 大 偏差 2.3%， 最 小 偏差 0.4%。 
计算 所 得 的 流量 值 与 实验 测量 结果 偏差 6.9%。 


一 一 1D network 
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图 8. 沿 程 压力 ( 气 源 0.5MPa) 
Fig 8.Pressure along the way(0.5MPa) 
如 表 4 为 气 源 压力 0.7MPa 时 ， 阀 门 不 同 开 度 工 
况 下 计算 流量 值 与 实验 结果 的 对 比 ， 开 度 02 工 况 下 
相对 偏差 最 大 ,为 10.47%,， 这 是 因为 小 开 度 下 流量 值 
小 ， 较 小 的 计算 偏差 就 会 造成 较 大 的 相对 偏差 ， 并 且 


0.3 0.966 1.074 11.15% 


0.4 1.209 1.325 9.57% 
0.5 1.419 1.526 7.56% 
0.6 1.519 1.655 8.96% 
0.7 1.605 1.720 7.17% 
0.8 1.626 1.738 0.88% 
0.9 1.652 1.768 7.05% 

1 1.660 1.788 7.68% 


如 表 6 为 气 源 压 力 0.35MPa 时 ,阀门 不 同 开 度 工 
况 下 计算 流量 值 与 实验 结果 的 对 比 ， 小 开 度 下 的 相对 
偏差 较 大 ， 开 度 0.3 工 况 条 件 下 最 大 为 11.8%。 综 述 
所 述 ， 网 络 计算 结果 与 实验 测量 结果 相对 偏差 绝 大 部 
分 工 况 下 小 于 10%， 少 数 工 况 小 于 13%， 验 证 了 空气 
网 络 模型 和 计算 方法 的 ;准确 合理 。 

表 6. 流量 对 比 (0. 35MPa) 
Table 6. Comparision of flow rate(0. 35MPa) 


实验 中 流量 测量 存在 一 定 的 误差 ， 管 路 中 较 少 量 气 流 
泄漏 至 环境 中 。 
表 4. 流量 对 比 (0. 7MPa) 
Table 4. Comparision of flow rate(0. 7MPa) 
Result of Result of 1D Relative 
exp(kg/s) network(kg/s) deviation 
0.1 0.555 0.566 2.13% 
0.2 0.900 0.995 10.47% 
0.3 1.334 1.360 1.96% 
0.4 1.617 1.704 5.39% 
0.5 1.886 1.956 3.74% 
0.6 1.977 2.141 8.30% 
0.7 2.095 2.227 5.94% 
0.8 2.185 2.270 3.89% 
0.9 2.264 2.306 1.84% 
1 2.234 2.330 4.32% 


如 表 5 为 气 源 压力 0.5MPa 时 ， 阀 门 不 同 开 度 工 
况 下 计算 流量 值 与 实验 结果 的 对 比 ， 在 小 开 度 下 的 相 
对 偏差 较 大 ， 开 度 为 0.2 下 最 大 为 12.24%。 
表 5. 流量 对 比 (0. 5MPa) 
Table 5. Comparision of flow rate(0. 5MPa) 


Result of Result of 1D Relative 
Value a 
exp(kg/s) network(kg/s) deviation 
0.1 0.410 0.430 4.86% 
0.2 0.626 0.702 12.24% 


Result of Result of 1D Relative 

exp(kg/s) network(kg/s) deviation 

0.1 0.293 0.326 11.37% 
0.2 0.447 0.500 11.75% 
0.3 0.671 0.750 11.80% 
0.4 0.905 0.924 2.17% 
0.5 1.034 1.039 0.49% 
0.6 1.088 1.108 1.82% 
0.7 1.161 1.176 1.30% 
0.8 1.210 1.202 -0.69% 
0.9 1.193 1.218 2.11% 
1 1.219 1.228 0.72% 

4 结论 


针对 某 实验 室 压缩 气体 供应 系统 ， 使 用 空气 网 络 
法 对 其 进行 研究 。 建 立 空气 网 络 模型 ， 通 过 实验 完善 
计算 方法 及 元 件 流动 特性 ， 进 行 空气 网 络 一 维 计算 。 
主要 结论 如 下 : 

1) 使 用 空气 网 络 法 研究 压缩 气体 供应 系统 ， 建 并 
专门 的 网 络 模型 。 空 气 网 络 中 有 阀门 、 局 部 损失 、 管 
道 。 提 出 新 型 数据 结构 应 用 于 空气 网 络 的 组 织 ， 图 形 
结构 -链表 -类 这 样 的 三 级 结构 ， 便 于 用 户 使 用 ， 并 且 
有 利于 节省 计算 机 内 存 。 

2) 阀 门 元 件 使 用 一 维 方法 计算 ， 通 过 实验 得 到 闪 
门 的 流动 特性 ， 流 量 系数 随 开 度 的 变化 曲线 。 并 将 计 
算 方 法 及 流量 系数 曲线 写 入 空气 网 络 程序 中 。 

3) 对 整个 空气 网 络 进行 一 维 计算 ， 计 算 结果 与 实 


Nl 


验 测 
压力 


确 可 


有 利 


量 结果 进行 对 比 。 气 源 压力 不 超过 0.7MPa 时 ， 

及 流量 相对 偏差 较 小 ， 从 而 验证 了 计算 方法 的 准 
言 。 针 对 压缩 空气 管 路 系统 的 一 维 计算 模型 研究 
于 较为 全 面 的 了 解 系 统 流动 特性 ， 实 验 中 对 阀门 


进行 快速 调控 。 


参考 文献 


[1] 


[2] 


[3] 


[5] 


[7] 


Yannick M. Secondary Air System Model For Integrated 
Thermomechanical Analysis Of A Jet Engine[R]. ASME 
Turbo Expo 2008: Power for Land, Sea and Air, 
GT2008-50078. 

Benra F K, Dohmen H J, Schneider O. Application of An 
Enhanced 1D Network Model To Calculate The Flow 
Properties of a Pre-swirl Secondary Air System[R]. ASME 
Turbo Expo 2008: Power for Land, Sea and Air, 
GT2008-50442. 

吴 丁 角 . 内 流 系统 的 网 络 计算 法 四 航空 学 报 ， 
1996,06:16-20. 

WU Dingyi. Network Technique of Intemal System[J]. Acta 


Aeronautica et Astronautica Sinica, 1996, 06: 16-20 (in 
Chinese). 

修 升 平 , 陶 智 等 . 非 稳 态 流体 网 络 方法 在 发 动机 空气 冷却 
系统 中 的 应 用 四 航空 动力 学 报 , 2009, 03, 494-498. 

HOU Shengping, Tao Zhi, et al. Application of Unsteady 


Fluid Network to Air Cooling System in Engine[J]. Journal of 
Aerospace Power, 2009, 03, 494-498 (in Chinese). 

郭 晓 杰 . 航空 发 动机 空气 系统 与 热 分 析 耦 合 方法 研究 [D]. 
上 海 交 通 大 学 , 2014: 10-17. 

GUO Xiaojie. 


Coupling Computational Approach of 
Secondary Air System and Heat Transfer in Aero-engine[D]. 
Shanghai: Shanghai Jiao Tong University 2014: 10-17 ( in 
Chinese) . 

吕 亚 国 , 刘 振 侠 等 . 航空 发 动机 内 流 空 气 系统 通 
个 设计 器 . 计算 机 仿真 , 2009,26(7):99-103. 

LU Yaguo, LIU Zhenxia. General Analysis Software Design 


分 析 软 


一 、 


of Aeroengine Intemal Air System[J]. Compuer Simulation 
2009,26(7): 99-103(in Chinese). 


胡 肖 肖 , 李 育 隆 等 . 基于 蒙特 卡 罗 法 的 航空 发 动机 稳 态 空 
气 系统 网 络 计算 加. 航空 发 动机 , 2014, Vol40, No.6, 28-32. 


HU Xiaoxiao, LI Yulong. Aeroengine Steady Air System 
Calculation Based on Monte Carlo Method[J]. Aeroengine, 
2014, Vol40, No.6, 28-32(in Chinese). 

潘 耘 峰 . 燃气 透 平 冷却 空气 系统 流体 网 络 法 研究 [D]. 
国 科 学 院 研究 生 院 , 2011: 53-56. 

PAN Yunfeng. Investigation on Fliud Network Method for 
The Analysis of Turbine Air Cooling System[D]. Beijing: 


UD 


9] 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


Graduate School of the Chinese Academy of Sciences, 2011: 
$53-56(in Chinese). 

陆 震 露 . 内 燃 机 一 维 流动 计算 方法 研究 及 性 能 仿真 软件 
设计 [D]. 上 海 交 通 大 学 , 2014: 15-25. 

LU Xiaoxiao，Computational Study on One-dimensional 


Fluid Flow and Desing of Performance Simulation Software 
for Internal Combustion Engine[D]. Shanghai: Shanghai Jiao 
Tong University 2014: 15-25(in Chinese). 

张 光 觅 , 许诺， 张 武平 . FLOWMASTER 在 暖 通 空调 中 的 
应 用 由. 制冷 与 空调 , 2006, No.03: 34-36. 

ZHANG Pengguang, Xu Nuo, Zhang Wu-ping. Application 
of FLOWMASTER in HVACI[J]. Refrigeration & Air- 
Condition, 2006, No. 03: 34-36(in Chinese). 


Meitner P L. Computer Code for Predicting Coolant Flow and 
Heat Transfer in Turbomachinery[R]. NA- SATP-2985, 1990. 
王 新 月 . 气体 动力 学 基础 [M]. 西安 : 西北 工业 大 学 出 版 社 . 
2006: 94-98, 204-209. 

WANG Xinyue. Fundamentals of Gas Dynamics[M]. Xian: 
Northwestern Polytechnical University Press. 2006: 94-98, 
204-209(in Chinese). 


Shapiro A H. The Dynamics and Thermodynamics of 
Compressible Fluid Flow Volume I[M]. New York: Ronald 
Press Co, 1954. 

吴 国 巾 . 调节 阀 使 用 与 
1999: 67-70. 

WU Guoxi. The Regulator Using and Maintenance[MI]. 


hi 


作 修 [TM]. 北京 : 化 学 工业 出 版 社 ， 


Beijing: Chemical industry press, 1999: 67-70(in Chinese). 
Ivor Horton， 苏 正 泉 ， 李 文 娟 ( 译 ). Visual C++ 2010 入 门 经 
(第 5 版 IIM]. 北 京 : 清华 大 学 出 版 社 , 2010: 459-539. 
Horton, 


Ivor SU Zhengquan, LI Wenjuan(translate). 
Beginning Visual C++ 2010(Fifth edition)[M]. Beijing: 


Tsinghua University Press, 2010: 459-539(in Chinese). 


